Influence of Grafting on Rootstock Rhizosphere Microbiome Assembly in Rosa sp. ‘Natal Brier’

Simple Summary Microorganisms are tightly associated with crops and can be pathogens or beneficials. Rose is the main ornamental crop worldwide. During production, rose varieties (over 100 in a single farm) are permanently grafted in a unique genotype as rootstock to improve plant performance. This work aimed to study the effect of grafting on root-associated microorganisms using next-generation DNA sequencing analysis. To this end, we have demonstrated that root-associated microorganisms of the rootstock (bacteria and fungi) will depend on the grafted genotype. In other words, a change in the variety will drive unique changes in the assembly of microorganisms at the root. This manuscript discusses differences in the bacterial and fungal communities when grafted and non-grafted and their potential impact on plant performance and agriculture. Abstract The root microbiome is vital in plant development and health and is highly influenced by crop cultural practices. Rose (Rosa sp.) is the most popular cut flower worldwide. Grafting in rose production is a standard practice to increase yield, improve flower quality, or reduce root-associated pests and diseases. ‘Natal Brier’ is a standard rootstock used in most commercial operations in Ecuador and Colombia, leading countries in producing and exporting ornamentals. It is known that the rose scion genotype affects root biomass and the root exudate profile of grafted plants. However, little is known about the influence of the rose scion genotype on the rhizosphere microbiome. We examined the influence of grafting and scion genotype on the rhizosphere microbiome of the rootstock ‘Natal Brier’. The microbiomes of the non-grafted rootstock and the rootstock grafted with two red rose cultivars were assessed using 16S rRNA and ITS sequencing. Grafting changed microbial community structure and function. Further, analysis of grafted plant samples revealed that the scion genotype highly influences the rootstock microbiome. Under the presented experimental conditions, the rootstock ‘Natal Brier’ core microbiome consisted of 16 bacterial and 40 fungal taxa. Our results highlight that the scion genotype influences root microbe’s recruitment, which might also influence the functionality of assembled microbiomes.


Introduction
Rose (Rosa sp.) is a woody shrub that belongs to the Rosaceae family and is the most economically important cut flower crop worldwide. In South America, cut-rose production is mainly located in Ecuador and Colombia. In 2019, the export value of Ecuadorian roses was nearly 800 million dollars, and their principal destinies were the USA, Russia, and the Netherlands [1]. Traditionally, Ecuadorian farms grow from 50 to 120 different rose cultivars in a single farm. Red cultivars usually occupy about 30% to 50% of farm Biology 2023, 12, 663 2 of 18 production land. Red roses are the most popular, and the cultivars Freedom TM by Rosen Tantau and Explorer TM by Interplant are the most popular red cultivars in Ecuadorian farms [2][3][4]. Grafting is a standard practice throughout commercial operations to increase yield, improve flower quality, or reduce root-associated diseases [5][6][7][8]. 'Natal Brier' (a hybrid of R. × multiflora and R. × damascena) has been the most used rootstock since 1990 [9,10].
Grafting is an ancient asexual propagation technique combining two plant genotypes to form a single composite plant. The composite plant has two sections, the scion (above-ground portion) and the rootstock (below-ground portion); each section has its genotype [11]. This technology widely propagates species that cannot be propagated by seed or cuttings. It also can improve growth, yield, and promote biotic/abiotic stress resistance [11]. Grafting has been reported to modify the metabolic profile of root exudates [12,13]. Root exudates play critical roles in plant nutrition, protection against pathogens, and microbiome assembly [14].
Thousands of microorganisms inhabit plants. The plant-associated microbial communities, their genomic material, and their interactions are known as the plant microbiome. Notably, the rhizosphere (soil zone closest to the root) is a hotspot for microbial activity influenced by root exudates [14]. The rhizosphere microbiome significantly modulates plant development and health [15] and can be influenced by phenological stage [16,17], biotic and abiotic stresses [18], or agriculture management practices [19,20]. A better understanding of microbiome dynamics in response to agricultural inputs and practices is vital to promoting the establishment of microbiomes that positively affect plant performance.
Few studies have looked at the effect of grafting on rhizosphere microbiome assembly. In Vitis vinifera, five rootstock genotypes grafted with the same scion genotype assemble different rhizosphere microbiomes; however, the observed rootstock genotype effect was farm dependent [21]. In Vitis vinifera cv. Lambrusco, potassium deficiency was associated with a rootstock genotype that assembles a microbiome depleted of potassium-solubilizing microorganisms [22]. Another study in grapevines showed that the scion cultivar highly influences the α-diversity and bacterial community structure at the rootstock rhizosphere and that the observed scion cultivar effect depends on the rootstock genotype evaluated [23]. Similarly, the apple scion genotype was reported to regulate root carbohydrate content and bacterial community structure and function [24]. In the case of herbaceous crops, it was reported that the rootstock genotype does not influence the leaf endophyte community assembly of tomato cv. 'Momotaro-Haruka' [25]. However, the structure of the rhizosphere microbiome of bottle gourd [Lagenaria siceraria (Molina) Standl.] was highly differentiated when grafted with a watermelon [Citrullus lanatus (Thunb.) Mansf.] scion. Interestingly, the rhizosphere microbiome shifted towards a structure similar to a non-grafted watermelon [26]. In the case of Rosa sp., prior studies have looked at endophyte communities across different genotypes [27] and root microbiota associated with rose replanting disease [28][29][30]. However, little is known about the influence of rose scion genotypes on the rootstock rhizosphere microbiome.
Due to the economic importance of grafting and the rootstock 'Natal Brier' for the Ecuadorian and Colombian rose industry, we aimed to understand the microbial composition and diversity at the rhizosphere of the rootstock 'Natal Brier' with and without grafting. To accomplish this objective, rhizosphere samples were collected from the Natal Brier rootstock with no graft and the Natal Brier rootstock grafted with two red rose cultivars (Explorer and Freedom). The 16S rRNA and ITS amplicon sequencing technology assessed bacterial and fungal community composition.

Sample Collection
In May 2020, the sampling was conducted in Arbusta 2 (latitude N 0 • 2 53.66 ; longitude O 78 • 11 30.85 ), a farm in Tabacundo, Ecuador (Figure 1a,b). A total of 12 bulk soil (100 gr each) and 12 root samples were collected. Four soil samples and four root samples were collected from two plant combinations, Natal Brier (NB)-Freedom and NB-Explorer, and the non-grafted rootstock (NB) (Figure 1). Explorer and Freedom are red bud color cultivars (Figure 1d,e). Sampled plants were randomly selected. Bulk soil samples were collected between the plants and 20 to 30 cm deep. For all soil samples, soil texture was classified as sandy loam with an average pH of 5.8 ± 0.3, 5 ± 0.2% of organic matter, and electrical conductivity of 2 mS ± 0.2. The roots were removed from the plant (rhizosphere samples) with a shovel at 10 to 15 cm depth, and the roots were gently shaken to remove loose soil. All soil samples were placed in polyethylene bags and stored in a portable cooler until transported to the laboratory.
In May 2020, the sampling was conducted in Arbusta 2 (latitude N 0°2′53.66″; longitude O 78°11′30.85″), a farm in Tabacundo, Ecuador (Figure 1a, b). A total of 12 bulk soil (100 gr each) and 12 root samples were collected. Four soil samples and four root samples were collected from two plant combinations, Natal Brier (NB)-Freedom and NB-Explorer, and the non-grafted rootstock (NB) (Figure 1). Explorer and Freedom are red bud color cultivars (Figure 1d,e). Sampled plants were randomly selected. Bulk soil samples were collected between the plants and 20 to 30 cm deep. For all soil samples, soil texture was classified as sandy loam with an average pH of 5.8 ± 0.3, 5 ± 0.2% of organic matter, and electrical conductivity of 2 mS ± 0.2. The roots were removed from the plant (rhizosphere samples) with a shovel at 10 to 15 cm depth, and the roots were gently shaken to remove loose soil. All soil samples were placed in polyethylene bags and stored in a portable cooler until transported to the laboratory.
The roots were placed in Falcon tubes filled with phosphate buffer solution (PBS) for rhizosphere samples and vortexed for 10 min to detach the rhizosphere area. Roots were removed and tubes centrifuged at 500 rpm for 10 min. The supernatant was discarded, and the samples were stored at −20 °C until further use [23].

DNA Extraction, Library Preparation, and Sequencing
Genomic DNA was extracted using the DNeasy Powersoil ® Isolation Kit following the manufacturer's instructions (Qiagen, Hilden, Germany). The DNA quality and concentration were assessed using a Qubit 4 (Invitrogen, Thermo Fisher Scientific, Waltham, MA, USA) and agarose gel electrophoresis. The roots were placed in Falcon tubes filled with phosphate buffer solution (PBS) for rhizosphere samples and vortexed for 10 min to detach the rhizosphere area. Roots were removed and tubes centrifuged at 500 rpm for 10 min. The supernatant was discarded, and the samples were stored at −20 • C until further use [23].

DNA Extraction, Library Preparation, and Sequencing
Genomic DNA was extracted using the DNeasy Powersoil ® Isolation Kit following the manufacturer's instructions (Qiagen, Hilden, Germany). The DNA quality and concentration were assessed using a Qubit 4 (Invitrogen, Thermo Fisher Scientific, Waltham, MA, USA) and agarose gel electrophoresis.
The ITS2 region was amplified using the primers ITS1-F (5 -CTTGGTCATTTAGAG-GAAGTAA-3 ) and ITS2 (5 -GCTGCGTTCTTCATCGATGC-3 ). DNA samples were diluted to concentrations of 3.5 ng µL −1 . Samples were run in triplicate. Each reaction contained 10 ng of DNA, 1 U of buffer, 0.3% bovine serum albumin, 2 mM of magnesium chloride, 200 µM of dNTPs, 300 nM of each primer, and 2 U of the DFS-Taq DNA polymerase (Bioron, Ludwigshafen, Germany). The PCR protocol was: 2 min at 94 • C, 25 cycles at 94 • C for 30 s, 55 • C for 30 s, 72 • C for 30 s, and 10 min at 72 • C to finish [31]. Then PCR products underwent enzymatic purification with exonucleases. A second triplicate amplification was performed with 3 µL of the purified PCR products and barcode-specific primers for the cycled samples under the previously mentioned conditions [31].
AMPure XP magnetic beads (Beckman Coulter, Brea, CA, USA) were used for bacterial and fungal amplicon purification and then quantified with the Qubit 2.0 fluorometer (Invitrogen, Carlsbad, CA, USA). Finally, the bacterial and fungal amplicons were bound in equal proportions and diluted to a concentration of 10 pM by sequencing, carried out on the Illumina MiSeq platform (Illumina, San Diego, CA, USA) using the 600-cycle V3 [31].

Data Analysis
16S rRNA and the ITS sequences were processed independently. MT-Toolbox software [32] was used to check sequence quality. Only sequence reads with 100% correct primer and a Q score over 20 were used. After the quality control, samples with less than 8000 16S rRNA reads or 17,000 ITS reads were discarded. The resulting sequences were clustered into amplicon sequence variants (ASVs) with the R package DADA2 version 1.8.1 [33]. The taxonomic assignment of each ASV was performed using the naïve Bayes k-mer method implemented in the DADA2 package using the Silva 132 database as a training reference. Fungal ITS forward sequence reads were processed using DADA2 [33] with default parameters. The taxonomic assignment of each ASV was performed using the naïve Bayes k-mer method implemented in the MOTHUR package [34] using the UNITE database [35] as a training reference.

Diversity Analysis, Community Composition, Community Function, and Statistical Analysis
16S rRNA and the ITS counts were rarified to 8000 reads per sample for bacteria and 17,000 reads per sample for fungi before further analysis with the 'ohchibi' package [36]. Alpha diversity indexes (Richness, Chao1, Shannon, Inverse Simpson, and Evenness) were calculated with the 'diversity' function from the vegan package version 2.5-3 [37]. Differences in the alpha diversity indexes between groups were analyzed with the Student's t-test (p value of <0.05) using R base functions [38].
Beta diversity was visualized through a principal coordinate analysis (PCoA) constructed using the Bray-Curtis dissimilarity index. Permutational analysis of variance (PERMANOVA) was performed using the adonis function of the vegan package version 2.5-3 [37].
The relative abundance of bacterial phyla and fungal classes were exemplified using the stacked bar representation encoded in the 'chibi.phylogram' function of the ohchibi package [36].
To compare the enrichment of specific ASV in non-grafted versus grafted rhizosphere communities, we used the 'DESeq2 package [39] to run the model abundance ∼ graft_status using the raw count table. An ASV was considered statistically significant if it had a false discovery rate (FDR)-adjusted p-value of <0.01.
The functionality of bacterial and fungal communities was analyzed using the 'Functional Annotation of Prokaryotic Taxa (FAPROTAX)' [40] and the 'FUNGuild' tools [41], respectively. These tools used the taxonomic assignation of each ASV and established its functionality based on specific databases. Significant differences between non-grafted and grafted were established with a Student's t-test (p-value of <0.05) using R base functions [38].
An ASV was considered part of the core microbiome if it was present in all samples (100% occupancy) with at least 5 reads. The sequence for the core ASVs were used to construct a maximum likelihood phylogenetic tree with the online tool 'RAxML BlackBox' [42]. We then used the Interactive Tree of Life (iTOL) web-based tool (https://itol.embl.de/) (accessed on 2 February 2021), to visualize the tree. 16S rRNA and ITS amplicon sequencing were performed on rhizosphere and bulk soil samples from the two plant combinations, Natal Brier (NB)-Freedom and NB-Explorer, and the non-grated rootstock (NB). 16S rRNA amplicon sequencing rendered 2107, 1514, and 1924 bacterial ASVs in NB, NB-Freedom, and NB-Explorer, respectively. Only 1588, 891, and 1302 bacterial ASVs were unique for NB, NB-Freedom, and NB-Explorer, respectively. ITS amplicon sequencing produced 676, 934, and 1028 fungal ASVs in NB, NB-Explorer, and NB-Freedom, respectively. Only 288, 483, and 542 fungal ASVs were unique for NB, NB-Explorer, and NB-Freedom, respectively. A total of 202 bacterial ASVs and 223 fungal ASVs were shared among all three rootstock systems (NB-Freedom, NB-Explorer, and NB).
Bulk soil and rhizosphere samples were separated using a principal component analysis (PCo2) (Figure 2). The rhizosphere and bulk soil assembled different bacterial communities ( Figure 2a). In contrast, the fungal communities in the rhizosphere and bulk soil samples clustered similarly (Figure 2c). There was no difference between the bacterial α-diversity of bulk soil and the rhizosphere (Table 1). However, we observed that bulk soil has higher fungal α-diversity than the rhizosphere (p-value < 0.05) ( Table 1). In soil and rhizosphere samples, the bacterial communities were dominated by the phyla Actinobacteria, Bacteroidetes, and Proteobacteria (Figure 2b), and the fungal communities were dominated by Ascomycota ( Figure S1) at the class level by Sordariomycetes, Dothideomycetes, and Mortierellomycetes ( Figure 2d).

Alpha Diversity, Taxonomic Composition, and Function of Microbial Communities Associated with the Grafted and Non-Grafted 'Natal Brier' Rootstock
The influence of the graft scion on microbial communities was evaluated in the rhizosphere samples. The rhizosphere bacterial and fungal communities in the non-grafted rootstock differed from those in the grafted rootstock. (Figure 3a,c). Bacterial α-diversity indexes were not statistically different but tended to be lower in the grafted rootstock ( Table 2). In contrast, fungal α-diversity was higher in the grafted rootstock (p-value < 0.05) ( Table 2). Ten bacterial phyla and twelve fungal classes dominated the rhizosphere microbiome of grafted and non-grafted plants. Proteobacteria was the most abundant bacterial phylum, followed by Bacteroidetes, Actinobacteria, and Acidobacteria. Sordariomycetes were the most abundant fungal class (~50%) (Figure 3b,d).   Differential abundance analysis was performed to identify microbial taxa significantly enriched in grafted and non-grafted plants (Figure 4a,b). The rhizospheres of non-grafted plants enriched seven bacterial ASVs assigned to the phyla Actinobacteria (ASV446), Alphaproteobacteria (ASV40), Bacteroidetes (ASV20 and ASV417), Gemmatimonadetes (ASV509), and Verrucomicrobia (ASV639 and ASV278) (Figure 4a    Differentially abundant bacterial and fungal taxa (ASVs) were assigned into functional categories using the FAPROTAX database [40] and the FUNGuild database [41], respectively. The abundance of bacteria involved in chemo heterotrophy, aerobic chemo heterotrophy, cellulolysis, and nitrate reduction was higher in non-grafted plants. On the other hand, the abundance of bacteria involved in nitrification, aerobic nitrite oxidation, intracellular parasites, and non-photosynthetic cyanobacteria was higher in grafted plants (Figure 5a). The abundance of fungal taxa classified as plant pathogens was higher in the non-grafted plants. In contrast, the fungal taxa in grafted plants were enriched into ectomycorrhizal and arbuscular mycorrhizal categories (Figure 5b).

Microbial Community Composition in the Rhizospheres of 'Natal Brier' Rootstock Grafted with Cultivars 'Freedom' and 'Explorer' and Their Core Microbiome
We observed an effect of the scion genotype on bacterial and fungal communities in the rootstock rhizosphere. Principal coordinate analysis showed that NB-Explorer rhizosphere samples separate from NB-Freedom samples along the first principal coordinate (PCo1) (Figure 6a,c). There were non-significant differences in bacterial and fungal α-diversity between plants grafted with Explorer and Freedom cultivars, but plants grafted with Explorer tended to have higher values on the α-diversity indexes evaluated (Table  3).
NB-Explorer and NB-Freedom plants shared 16

Microbial Community Composition in the Rhizospheres of 'Natal Brier' Rootstock Grafted with Cultivars 'Freedom' and 'Explorer' and Their Core Microbiome
We observed an effect of the scion genotype on bacterial and fungal communities in the rootstock rhizosphere. Principal coordinate analysis showed that NB-Explorer rhizosphere samples separate from NB-Freedom samples along the first principal coordinate (PCo1) (Figure 6a,c). There were non-significant differences in bacterial and fungal α-diversity between plants grafted with Explorer and Freedom cultivars, but plants grafted with Explorer tended to have higher values on the α-diversity indexes evaluated (Table 3).

Discussion
Rose grafting is a propagation technique widely used to improve plant performance and flower quality [7,8,43]. The superior performance of grafted rose plants depends on rootstock selection, growing media, geographic region, environmental conditions, and culturing practices [7,[44][45][46]. The rootstock genotype selection influences the reblooming time, scion nutritional status, and dry matter partitioning [47,48]. Several reports studying rose scion-rootstock interaction are available; however, to our knowledge, none have focused on the rootstock microbiome and the influences of scion genotype. We described the rhizosphere microbiome of the non-grafted rootstock 'Natal Brier' and how it changes upon grafting with red cultivars 'Explorer' and 'Freedom. The plant and its associated microbiota are thought to function as an entire entity. The plant-associated microbiota serves the plant by providing additional processes that facilitate or improve plant development, health, and adaptation [49]. Understanding the effects of grafting on the rootstock microbiome is necessary to unveil the mechanism behind the superior performance of grafted plants. The bacterial and fungal communities of the rootstock' Natal Brier' are dominated by the phyla Proteobacteria and Sordariomycetes, respectively (Figure 3b,d). Similarly, the bacterial microbiome of Rosa corymbifera 'Laxa', a rootstock sensitive to rose replanting disease (RRD) used for garden rose propagation, is dominated by the phylum Proteobacteria. However, the soil source determines the dominant phyla of the fungal microbiome of R. corymbifera 'Laxa' [28]. Thus, the soil source is an important factor to consider for further understanding of rootstock microbiome assembly.

Grafting Changes the Rhizosphere Microbiome of Rose sp. 'Natal Brier'
Root exudates play vital roles in plant nutrition, protection against pathogens, and microbiome assembly [14]. Upon soil colonization, root compartments are occupied by soilborne microorganisms. Primary and secondary metabolites regulate plant-microbe interactions at the root zone into the rhizosphere [15]. Shifts in the metabolic profile of root exudates are associated with plant developmental stage, health status, or stress responses [50]. In addition, the changes in root exudate composition are accompanied by changes in microbiome structure and diversity [17,51]. In R. corymbifera 'Laxa', the rhizosphere microbiome and root metabolite content change in response to rose replant disease (RRD) [28]. Grafting induces changes in microbiome structure, microbial activity, and root exudate composition of the rootstock [12,13,26,52]. Consistent with previous observations in other plant species, grafting also changed the rhizosphere microbiome composition of the rootstock 'Natal Brier' (Figure 3). Other reports have shown that bacterial diversity in the rhizospheres of grafted plants is higher than that of non-grafted rootstock rhizospheres [53]. However, we observed that grafting only increased the diversity of the fungal microbiome but not bacterial diversity ( Table 2).
Bacterial functional prediction in grafted plants included nitrification, aerobic nitrite oxidation, intracellular parasites, and non-photosynthetic cyanobacteria. Accordingly, the genus Nitrosospira and Nitrospira were enriched. Nitrosospira converts ammonia (NH 3 ) to nitrite (NO 2− ) and Nitrospira converts nitrite to nitrate (NO 3 ) during the nitrification process [54]. The rhizospheres of grafted plants was enriched with nitrogen-fixing taxa (Frankiales and Rhizobiales) [55,56]. Other beneficial bacterial taxa included Chitinophagaceae and Saccharimonadales [57,58]. The increase in the abundance of nitrogen-fixing bacteria in grafted plants suggests a mechanism for superior performance in grafted plants. Likewise, bacterial taxa involved in nitrogen cycling can significantly impact nitrogen use efficiency.
Conversely, the functions assigned to bacteria enriched in the non-grafted rootstock were chemoheterotroph, aerobic chemoheterotroph, cellulolysis, and nitrate reduction. For example, Jatrophihabitans soli KIS75-12T was reported to hydrolyze urea into ammonia [59]. In addition, Gemmatimonadaceae abundance correlates with N fertilization [60]. Other bacterial taxa reported as beneficial are Sphingomonadaceae and Verrucomicrobia [61]. Species of genus Niastella have been isolated from soil, and Niastella gongjuensis sp. nov was reported to reduce nitrate [62].
The fungal taxa enriched in the rhizospheres of the non-grafted plants were assigned to the plant pathogen category. These taxa were classified as Alternaria spp., Blastobotrys spp., Chaetomiaceae spp., Cladosporium spp., Disculoides spp., Fusarium spp., and Plectosphaerella spp. Among them, only Cladosporium spp., Alternaria spp. and Fusarium spp have been previously reported to cause disease in roses. Cladosporium cladosporioides causes petal discoloration and necrosis [63], Alternaria spp. causes bud blight and leaf black spot [64], and Fusarium oxysporum causes dry rot diseases in rose seedlings [65].

Scion Genotype Determines the Changes in the Rhizosphere Microbiome of Rootstock Rose sp. 'Natal Brier'
The performance of rose cultivars grown on their own roots is inferior to that of plants grafted onto a rootstock. The superior performance of grafted roses results from the additive effect of the scion and rootstock genotypes. Thus, the superiority of a grafted rose onto a specific rootstock depends on the scion vigor [8]. Other studies showed that the scion genotype could modify root biomass, transcriptome, and root exudate composition [87][88][89]. Likewise, the scion genotype can mediate microbiome recruitment [52]. To gain more insight into the scion genotype effect on the rhizosphere microbiome, we analyzed the microbiome of the red cultivars 'Freedom' and 'Explorer' grown onto the rootstock 'Natal Brier'. Scion genotype determined the rhizosphere microbiome composition ( Figure 6). Although no significant differences were found, both fungal and bacterial alpha diversity of the cultivar 'Explorer' tended to be higher than of the cultivar 'Freedom'. The soil has been highlighted as a reservoir of microbial diversity [90], and the root exudates drive microbiome recruitment. As previously discussed, the scion genotype is an essential factor in the composition of the root exudates and rhizosphere microbiome. Accordingly, the rhizosphere microbiome of the grafted rootstock 'Natal Brier' showed taxa unique for each rootstock-scion combination, and the group of taxa consistently present is considered as a core microbiome (Figure 7). The number of unique microbial and fungal taxa was higher in the rootstock grafted with the cultivar 'Explorer'. The bacterial core microbiome consists of taxa belonging to the phyla Proteobacteria, Acidobacteria, Actinobacteria, Firmicutes, and Patescibacteria. The fungal core microbiome consists of taxa belonging to the phyla Sordariomycetes, Eurotiomycetes, Dothideomycetes, Ascomycota_unclassified, Orbiliomycetes, Tremellomycetes, Tremellomycetes, Archaeorhizomycetes, Mortierellomycetes, and Malasseziomycetes. Further investigation is needed to reveal the potential contributions of the microbiome to the yield and vase life of these red cut-rose cultivars.

Conclusions
Grafting is a propagation technique that is widely used to improve productivity and flower quality in cut-rose production. The superior vigor of grafted rose cultivars results from the additive effect of the rootstock and scion. In grafted plants, the scion may mediate microbiome assembly in the rootstock [52]. Here, we described the effect of grafting on the rhizosphere microbiome of the rootstock 'Natal Brier'. Our data suggest that the rhizosphere microbiome of the rootstock changes upon grafting. Grafting leads to the enrichment of plant-growth-promoting taxa in the rhizosphere. Changes in microbial community composition also modulated the microbiome functions. A wider range of microbiome functions potentially contributes to the superior performance of grafted plants. According to our data, the scion cultivar modulates the recruitment of microbiomes with beneficial effects. Thus, assessment of a wider diversity of scions will contribute to discovering microbiomes with superior beneficial functions. Understanding the potential functions of the plant microbiome opens an avenue to exploit the full potential of crops by leveraging the microbiome through adapting current breeding strategies, agronomic practices, and microbial products. The rose rootstock' Natal Brier' is a suitable biological model for basic and applied research to study the effects of grafting on the microbiome due to the vast number of scion genotypes available and its economic importance for countries cultivating cut roses [10].  Data Availability Statement: Amplicon sequencing data are available at the NCBI Sequence Read Archive (project PRJNA893836). Count tables and relevant data files can be found at GitHub https: //github.com/darioxr/grafting_rose_microbiome (accessed on 10 December 2022).